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Abstract 

It is shown that quasi-Langevin method can be used for the calculation of the 
shot noise in correlated single-electron tunneling. We generalize the existing 
Fokker-Plank-type approach and show its equivalence to quasi-Langevin ap- 
proach. The advantage of the quasi-Langevin method is a natural possibility 

to describe simultaneously the high ("quantum") frequency range. 
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Correlated single-electron tunneling |l|] remains an attractive topic during last decade. 
Because in the systems of small-capacitance high-resistance tunnel junctions electrons tunnel 
almost as the classical particles, the most of experiments are well explained within the 
framework of "orthodox" theory [|I| of single-electronics which is based on the classical 
master equation. 

The theory of the shot noise in single-electron transistor p| developed in Refs. is also 
based on the master equation. The shot noise appears due to the randomness of the tunneling 
events. Despite the classical description of the system, the current in this theory should be 
treated as a kind of an operator because the current is caused by tunneling events which 
change the charge state of the system. This theory has been studied in detail in a number of 
subsequent publications - see, e.g., Refs. 0-0. The first experimental confirmation of the 



theory has been obtained [|T2| in 1995, and other experiments are definitely coming because 
several groups has succeeded in the fabrication of the single-electron transistor operating at 
relatively high frequencies (hundreds of kHz) at which 1/f noise should become less than the 
shot noise. 

Among the most noticeable developments of the shot noise theory in single-electronics 
after 1991, let us mention the following. It has been shown that the noise due to the 
processes of Andreev reflection which transfer two electrons per tunneling event can be 
described by the same formalism with e substituted by 2e. The shot noise in single-electron 
transistor with large level discreteness has been studied |13[ . The shot noise theory has been 
applied to single-electron systems other than single-electron transistor [^4|,|15|. Besides the 
noise in the "orthodox" frequency range uj ~ W/ e^R (where is a typical energy and R is 
a typical resistance), the noise in the "quantum" frequency range u ~ W/h has been studied 



and the matching of the high and low frequency limits of two ranges has been proven P 
however, the approach unifying both frequency ranges in one formalism has not been found. 

The existing formalism of the shot noise in single-electron tunneling is of the Fokker- 
Plank type and is based on the deterministic master equation. However, generally more 
popular method in the study of the noise is the Langevin approach in which the random 



term is introduced into the evolution equation. Three decades ago the Langevin approach 
was apphed [|I6| to study the fluctuations in the nonequilibrium electron gas. In the present 
letter we show that a similar method can be used for the single-electron tunneling. (We call 
it quasi-Langevin because the regular Langevin method cannot be applied directly.) It is 
equivalent to existing formalism in the "orthodox" frequency range, but it naturally allows 
also the calculations in the "quantum" frequency range. 

Let us start with the generalization of the existing Fokker-Plank type method to an 
arbitrary system consisting of voltage sources, capacitances and tunnel junctions with suf- 
ficiently large resistances [Rj ^ Rk = h/e^ ~ 2QkVL). For simplicity we do not consider 
Ohmic resistances. In this case the dynamics is governed by the matrix master equation 

a = Ta (1) 

where the vector (T„ is the probability to find the system in the charge state n = {rii, ...ni} 
(which is characterized by the numbers Ui of excess electrons in each of L internal nodes of 
the system) and 

k j 

where L^^^ are the corresponding tunneling rates and the summation over the junction 
number j is necessary when an electron can tunnel to (from) an internal node from (to) 
different external electrodes. 

To find the mutual spectral density for two processes X{t) and Y{t) we can calculate first 
the correlation function Kxy{'^) = {X{t ~^ '^)Y{t)) " {-^)0^) (brackets denote the averaging 
over time) and then take the Fourier transform Sxri^^) = 2 KxriT) exp(iujT) dr . If 
both X and Y are functions of the charge state n (for example, potential of a node) then 
the correlation function is given by the simple expression 

KxY{r) = e{T) X{m) a(r, m\n) Y{n) 
+ e{-T) Y: Y{m) a(m, -r|n) X(n) - (X) (F) , (3) 



where o"(r, m\n) is the retarded Green's function of Eq. ([^) being the probabihty to find the 
system in the state m at t = r > if at t = it was in the state n, {X) = I]n^(^)'^n*! 
0"^* is the stationary distribution, Fcr** = 0, J2n'^n — 1- (Notice that X and Y are classical 
variables, and their commutator is zero.) 

However, if X and/or Y represent the current through a tunnel junction or in an external 
lead, the Eq. (§) should be modified. For example, if X{t) is the current contribution 
corresponding to tunneling events T^^^^ while Y{t) corresponds to r:^/^„/, then (similar to 
Refs. 

KxY{r)/eiei = ^(r) rf„^„ a(r, n\m') 

+ ei-r) r£,^„,a(n', -r|m) F„ 

- r4^„<* ^C^nK' + Smm'5nn'5jj'5 (t) . (4) 

Here the last term is responsible for the high-frequency limit. The effective charges and 
€■£ are determined by the direction of electron tunneling, e-?^ = —e-L, and by the circuit 
capacitances (so that = e only if the current through junction j is measured). Any 
current-current correlation function can be written as a sum of Kxy{t) given by Eq. (^ 
over all possible transitions between charge states (such a sum is a counterpart of Eq. (^ 
in which the sum is written explicitly). The processes of cotunneling [|T7| in the simple 



master equation approximation []TB[ and fluctuations due to Andreev reflection |^ can be 
accommodated in this technique using the appropriate values of (for cotunneling the index 
j should obviously represent the set of junction numbers). 

For the correlation functions when X is a current and F is a function of the charge 
state (or vice versa), the recipe is the "combination" of Eqs. (^ and (^ while the term 
proportional to 6{t) is absent. 

The expressions for spectral densities directly follow from Eqs. (^) and (H) because the 
Fourier transformation affects only the evolution operator a{T,m\n), and the corresponding 
Green's function in the frequency representation is simply obtained from Eq. ([l|) 

a(iu,m\n) = (—iujl — T)~^ , (5) 

L J mn 
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where 1 is the unity matrix. For example, Eq. (|) leads to the following spectral density 



5xyM/e4e4 = 2r: 

^3 



+ 2P , 

' m'i—n 



-rl St 
nm' ^'^ri"^n' 



nm i — ri^ fn 



n'm 



+ 2 6nn'Smm'Sjj' T^rn^^ ■ (6) 

This method allows to calculate all spectral densities within the framework of "or- 
thodox" theory, and at least for the single-electron transistor the numerical procedure is 
rather trivial because matrix T is three-diagonal and the calculation of Eq. (Bf) is straitfor- 



ward [|T^]. There are some additional simplifications for — > 0, however, in this paper we 
concentrate on the finite frequencies only. Notice that in the special case when there are 
only two charge states involved, the theory is equivalent to the noise calculations in resonant 



tunneling diode pO| 



Now let us develop the Langevin-type approach. It is obvious that for the discrete 
random dynamics of the charge states of a single-electron circuit it is impossible to introduce 
the random term in Eq. (|I]) in a reasonable way. However, let us imagine the ensemble 
of M (M ^ 1) independent similar circuits, and let us average all magnitudes over this 
ensemble. Then the average (over time) currents and voltages will not change (due to 
ergodicity), but the spectral densities of fluctuations (second order magnitudes) will decrease 
M times. Hence, to calculate spectral densities of the initial system, we can take the leading 
(~ M^^) order of the spectral density of magnitudes averaged over the ensemble, in the 
thermodynamic limit M — ^ oo. In contrast to the single system, the dynamics of the 
ensemble at M oo can be easily described by the Langevin approach because there are 
no longer step-like dependencies on time. 

Now Mcr„(t) which is the number of systems in the ensemble having the charge state 
n, can be considered as non- integer (because of M ^ 1) and its random dynamics can be 
described with the help of fictitious random term ^(t). In the "orthodox" framework (when 
we consider electron jumps as instantaneous events), ^(t) is 5-correlated (white), and its 
amplitude is defined by the Poissonian nature of the tunneling process. The recipe is the 



following |T^: for each average flux MTi^^^a^ (number of transitions per second) in the 
space of charge states, we should add in the master equation the random flux ^l^^^it) with 
the corresponding "seed" spectral density given by usual Schottky formula 



n 



S.J .j' (^) 



2M ^5rr,.rn'S.r,r,'Sii'Ti„. 



(7) 
(8) 
(9) 



Notice that when there are fluxes in opposite directions {m ^ n and n <— m), we should 
apply C, {t) for each direction, so that the random flux does not vanish even if the net average 
flux is zero. 

Because of the linearity of Eqs. (0)-® the final spectral densities of the averaged (over 
M) magnitudes are obviously proportional to 1/M. Hence, rescaling to the single system can 
be done formally assuming M = 1 in Eqs. (0)-(|^). So, instead of keeping M and rescaling 
at the final stage, we will use M = 1 in all equation below. We call this simple trick of 
considering first the large ensemble, writing the Langevin equation for it, and then returning 
to a single system, a quasi-Langevin approach. 

Using the standard procedure we find the Fourier transform 



^n(^) 



(10) 



Then for the occupation-occupation spectral density 



n' 



(11) 



Using Eq. (P) after simple algebra we obtain 



'? — 9 
+ 2 



a, 



St 



cr, 



St 



(12) 
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which coincides with the result of Fokker-Plank approach (Fourier transform of Eq. 
without X and Y factors). 

The technique is similar for the current-current fluctuations. The case of Eqs. and 
(0) corresponds to currents 



X{t) 



Y{t) = ~ei 



m'*—n 



(13) 



and the straitforward (though rather lengthy) calculations using Eqs. (|) and (|l0|) lead to 
Eq. (D. 

The final expression for the current-occupation spectral density is 



+2 



'J 



[iuol 



St 



km 



nk 



m-i—n^n ' 



(14) 



and it also coincides with the corresponding expression obtained in the Fokker-Plank tech- 
nique. 

Thus, we have proven that the Fokker-Plank method is equivalent to the quasi-Langevin 
method within the "orthodox" framework. However, in contrast to the former approach, the 
quasi-Langevin method easily allows generalization for the fluctuations in the "quantum" 
frequency range. 

Let us remind that in "orthodox" theory (from now on for simplicity we speak only 
about purely single-electron tunneling and do not consider cotunneling, Andreev reflection, 
etc.) 

loiW/e) 



W = eV- e72C, 



e(l-exp(-iy/T))' ^' ^ ^^^^ 

where /o('^^) is the "seed" /-F curve of the junction (in the linear case Io{v) = v/R), W is 
the energy gain due to tunneling, V is the voltage across the junction before the tunneling, 
and Ceff is the effective junction capacitance (which also accounts for the environment). 
The generalization of quasi-Langevin method is the substitution of Eq. (|^) by more general 
equation ||14| 



St 



IoAW^^n/e±nu;/e) 



T- (16) 

e[l-exp[{Wl^n±huj)/T)\ 

This equation can be considered as a generalization of the fluctuation-dissipation theorem 



and equations of Ref. |^T| for the case of single-electron tunneling. It can be obtained 
within the standard tunneling hamiltonian technique averaging the quantum current-current 
correlator and then taking the Fourier transform. 

Using this generalization one can see that at high frequencies, u ^ T, the occupation- 
occupation and occupation-current spectral densities vanish, while for the current-current 
spectral density Eq. @ transforms into 

Sxy{uj) = eie^S , f (cj), (17) 

(given by Eq. (|T6|) ) because the first terms of Eq. (|^) are too slow to give a cotribution. This 



result coincides with the result of Ref. [0. The advantage of the quasi-Langevin approach 
is the possibility to obtain spectral densities in the "orthodox" and "quantum" frequency 
ranges using the same formalism while in Ref. [|1^ they are necessarily treated on different 
footing. 

In the "quantum" frequency range the current spectral density does not correspond 
directly to the available power because of the contribution from zero-point oscillations. The 
spectral density of the current calculated above can be considered as the power (within the 
unit bandwidth) going from the system to a small external resistance r (divided by r and 
the coupling factor a). To obtain the available power, we should subtract the power flow in 
the opposite direction which is the product of the voltage spectral density of the resistance 
2%ujr coi\i{huj / 2Tr) (the "receiver" temperature can differ from T), the factor a, and 
the active conductance of the fluctuation source, which in the case corresponding to Eqs. 
(|1|)-([13) is given by [|14 



ReG{u) = ({^^f/2hu) [f + - f-] (18) 
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(this expression obviously corresponds to the lowest order of photon-assisted processes). 
Hence, the contribution to be subtracted from 5*^ is equal to 2hLj coth.(hu/2Tr)ReG(uj). 
Traditionally this contribution is called zero-point for = 0. However, for finite tempera- 
ture T a more natural choice of the "receiver" temperature is Tj. = T. Notice that to get the 
total conductance ReGt{uj), Eq. ( [T^ ) should be summed over all kinds of tunneling events. 

The solid line 1 in Fig. |l] shows the numerical result for the spectral density Sii{uj) of 
the current in the external lead of the single-electron transistor consisting of two similar 
tunnel junctions with capacitances Ci = C2 = Cs/2 and resistances Ri = R2 = -Rs/2 at 
Qo = 0.3e and T = O.OSe^/Cs. The frequency dependence in the "orthodox" frequency 
range is important at / = uj/2tt ~ I/e while the "quantum" frequency dependence occurs 
at / ~ e^/CY.h. Because the junction resistances are chosen sufficiently large, i?s = IOORk, 
two frequency ranges are far away from each other (the ratio of typical frequencies is on the 
order of R/Rk)- The line 2 shows the spectral density corresponding to the power available 
at zero-temperature receiver (zero-point contribution is subtracted, = 0) while the line 3 
shows the available power for = T. Notice that the curve 3 differs considerably from curve 
1 even in the "orthodox" frequency range because of thermodynamical reasons. The dashed 
lines 4 and 5 show the results for T = (two lower curves obviously coincide). Notice 
the cusps at hf = 0.05e^/Cs which correspond to the extra energy above the Coulomb 
blockade and the abrupt vanishing of the available power above the frequency / = Wmax/h 
(/ = OAbe^ /Cj:,h in this case) where Wmax is the maximal energy gain among all possible 
tunneling events. 

In the case Rj S> Rk the "orthodox" and "quantum" frequency ranges are far apart 
from each other because F ~ W/eR <^ W/h. The fact that the quasi-Langevin formalism 
describes the fluctuations for the whole frequency axis, suggests the possibility to use this 
method also when two frequency ranges are close to each other. However, this is impossible 
because when Rj ~ Rk the master equation approach fails due to strong cotunneling pro- 
cesses |]T^. The correct description of the noise in this case (as well as just calculation of 
the average currents) is still an open question despite an important progress for the average 

9 



currents 122 



In conclusion, we developed the quasi-Langevin method for the calculation of the fluctu- 
ations in single-electron tunneling. Its equivalence to the previous Fokker-Plank approach 
in the "orthodox" frequency range is proven. The advantage of the quasi-Langevin method 
is the natural generalization for the "quantum" frequency range. 
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S. Stepanov, and J. S. Tsai for the valuable discussions. The work was supported in parts 
by US AFOSR, Russian Fund for Basic Research, and Russian Program on Prospective 
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FIGURES 

FIG. 1. Curve 1 shows the spectral density of the current (in the external lead) through 
the symmetric single-electron transistor at T = 0.03e^/Cs as a function of frequency / = uj/2t:. 
The relatively large junction resistances provide cosiderable difference between "orthodox" and 
"quantum" frequency scales. The curve 2 shows the spectral density corresponding to the available 
power with zero-point cotribution subtracted (T^ = 0) while for the curve 3 the "receiver" is 
assumed to be at the same temperature (Tj. = T). Dashed curves 4 and 5 demonstrate the results 
for T = 0. The average current / = 0.106e/i?sCE for T = O.OSe^/Cs while / = 0.09e/i?sCE for 
T = 0. 
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